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Self Organized Criticality in an one dimensional 
magnetized grid. Application to GRB X-ray 

afterglows 
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Abstract 

A simplified one dimensional grid is used to model the evolution of 
magnetized plasma flow. We implement diffusion laws similar to those 
so-far used to model magnetic reconnection with Cellular Automata. As 
a novelty, we also explicitly superimpose a background flow. The aim is 
to numerically investigate the possibility that Self-Organized Criticality 
appears in a one dimensional magnetized flow. The cellular automaton’s 
cells store information about the parameter relevant to the evolution of 
the system being modelled. Under the assumption that this parameter 
stands for the magnetic field, the magnetic energy released by one volume 
during one individual relaxation event is also computed. Our results show 
that indeed in this system Self-Organized Criticality is established. The 
possible applications of this model to the study of the X-ray afterglows of 
GRBs is also briefly considered. 

Keywords: radiation mechanisms: non - thermal, gamma - ray burst: gen¬ 
eral, X - rays: bursts. 


1 Introduction 

GRBs are explosive cosmic gamma-ray emissions, with typical energy fluxes of 
the order of 10“® to 5 x 10“^ J m“^, and durations that range from 10“^ to 10^ s 
[32] . Since their distribution is isotropic, they are believed to have a cosmological 
origin, being located at extra-galactic distances [32]. The observable effects of 
GRBs are produced by the dissipation of kinetic energy from a relativistically 
expanding plasma, though the underlying progenitor model is, as yet, unknown. 
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Reviews and discussion of new models may be found in, e.g., n m 113 SOI 1 

EHiiTiisniEiiiini. 

The detection of X-ray flares in nearly half of the observed GRB afterglows 
has added a new mystery to the GRB puzzle. The present astrophysical 
observations suggest that the central engines of bursts, after the gamma-ray 
emission has ended, still have long periods of activity, during which energetic 
explosions eject relativistic materials, leading to late-time X-ray emission |1H] - 
The GRB observations show that strong magnetic fields play a key role during 
the acceleration of the outflow and the prompt emission [23] , while any remain¬ 
ing magnetization at large distance can affect the interactions of the flow with 
the external medium. 

A statistical analysis of the X-ray flares of GRBs with known redshifts was 
performed in [48]. It was shown that X-ray flares and Solar flares share three 
statistical properties: a) power-law energy frequency distributions, b) power-law 
duration-time frequency distributions, and c) power-law waiting time distribu¬ 
tions, respectively. All these distributions can be well understood within the 
physical framework of a magnetic reconnection-driven Self-Organized Criticality 
system. These statistical similarities, together with the fact that Solar flares are 
triggered by a magnetic reconnection process taking place in the atmosphere of 
the Sun, suggest that X-ray flares originate from magnetic reconnection-driven 
events possibly involved in ultra-strongly magnetized millisecond pulsars |48j . 

In [TU] and [17] it was already proposed that GRB X-ray flares may be 
powered by magnetic reconnection events. Magnetic reconnection may occur 
near the photosphere if the outflow develops an alternating field structure due to 
e.g. magnetic instabilities or a misalignment between the magnetic and rotation 
axes. The magnetic dissipation model is favoured by observations m- 

13 proposed that the numerous spatial extended systems, which exhibit a 
number of properties which may be shortly characterized as flicker noise for the 
temporal evolution and self-similar (fractal) behavior for the spatial evolution 
can be organized under the same principle of Nature, namely as systems in the 
state of Self Organized Criticality (SOG) [55] - 

There is a lot of literature covering the aspects of SOG in astrophysics (see, 
e.g., 0 and references therein), ranging from simple but effective numerical 
models EnmEn], to sophisticated analytical models as in [52| and m- The 
possible areas of physics in which these models prove to reproduce observational 
data is very wide. This is a consequence of the fact that SOG may be an 
underlying principle in Nature, as the initial intuition of |3 stated. Examples 
in astrophysics include organization of plasma in accretion discs [28] to explain 
observational data such as those discussed in (33], axion clouds around black 
hole [30|, solar flares [giiiiiniEii. 

[39] have analyzed GRB X-Ray data, and they did find that all observational 
results can be fitted into a framework based on a stochastic pulse avalanche 
model running in a near-critical regime. Their reasoning is reproduced here: 
1. All GRBs can be described as different random realizations of the same 
simply organized stochastic process within narrow ranges of the parameters of 
the process. 2. The stochastic process should be scale invariant in time. 3. 
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The stochastic process works near its critical regime. This would explain the 
large morphological diversity of GRBs. They propose that the flaring events 
are based on magnetic reconnection, but their simulations do not include evo¬ 
lution equations for magnetic helds. While this type of statistics of data can be 
explained by turbulence, the analysis performed by |48j is clearly in favour of a 
Self Organized Criticality model. 

It is the goal of the present paper to investigate from the point of view of SOC 
one dimensional magnetized systems undergoing slow driving. More specifically, 
this study aims at describing a specific behaviour, i.e., that of magnetic field on 
a space scale that is global. Thus we do not attempt to solve the MHD equations 
nor do we even take all of them into account. We focus on the magnetic field as 
a being a relevant parameter, since it is a field that couples to the radiation field 
and since radiation is our main source of information on astrophysical sources. 
Also, the magnetic field is a parameter that exhibits a transition in its behaviour 
and this transition can be explained based on the Laplacian of the magnetic 
field reaching a critical value. The evolution equation for this parameter is the 
magnetic induction equation, which describes the time behaviour of a magnetic 
field in a resistive plasma medium in which the MHD approximation holds |d4] , 
with a background flow. 

Magnetic reconnection was investigated independently of the idea of SOC 
as a natural phenomenon occurring in all branches of astrophysics [sum]- It 
can be viewed as a change in the magnetic field topology when the value of 
a certain parameter locally reaches a threshold value. The change in topology 
occurs fundamentally because the configuration is no longer of minimum energy. 
A part or all of the excess energy is lost to radiation. Extension of the idea 
of magnetic reconnection to SOC is nothing more but moving from the local 
scale of the reconnection site to the entire system. So conceptually SOC in a 
magnetized environment is expected and it is thought to be the answer to many 
unanswered questions in astrophysics [4S|3]. This expectation was confirmed 
by observational data of solar flares (a review of theory and data is presented 
in, e.g., [7]). 

The present paper is organized as follows. In Sectionj^we introduce the mag¬ 
netic induction equation for plasma flows, and we present its one-dimensional 
version. The discretization of the induction equation as well as the simulation 
procedure is presented in Section The results of our simulations for one di¬ 
mensional magnetic flows and the establishment of the SOC are presented in 
Section The applications of the model for the study of CRB afterglows are 
briefly presented in SectionWe discuss and conclude our results in Sections]^ 
and[3 

2 Self organized criticality in magnetized flows 

The aim of this section is to investigate whether or not SOC is possible in a one 
dimensional magnetic field set-up. 
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2.1 Preliminaries: The theoretical model 

The set of MagnetoHydrodynamic equations which (under some assumptions 
regarding the ratio of microscopic and macroscopic time-and-space scales) ana¬ 
lytically describe general magnetized flows are inhomogeneous nonlinear partial- 
differential equations, which are sometimes anisotropic. Numerical techniques 
to solve this set of equations have evolved, but one still cannot solve them in 
full generality not even locally. 

Cellular Automata (CA) simulations are a valid alternative to the full nu¬ 
merical treatment when only qualitative and quantitative answers of a specific 
kind are sought for. In such a setup the velocity profile has not been consistently 
included, to our knowledge. In fact, most of the times it has been considered 
nonexistent and the convective part was replaced by a stochastic loading phase. 
The following treatment aims at a middle ground between ignoring the flow 
completely and the exact numerical solution of the equation (not accessible to 
a CA of this type). 

The skeleton of our approach to calculating the stress, number of events 
and released energy has been used in the literature, see e.g. EH; the novelty in 
this work is explicitly taking into account the background flow characteristic to 
GRBs. 

Thus we consider the induction equation 

^ = Vx (uxr) (1) 

and define the control parameter as 

G = —(2) 
nn 

where B is the magnetic field, v is the flow velocity, rj is the constant plasma 
resistivity and the notation nn denotes the number of nearest neighbours in the 
conhguration. 

The set-up consists of a flux tube aligned with the Oz axis of a system of 
coordinates fixed at the footpoint of the tube. A magnetic held is present, B — 
{B{z),0,0) and we allow for a background how v = (0,0, w(z)). We emphasize 
that the setup is thus two dimensional. In this setup, with nn = 2, the magnetic 
induction equation and the control parameter become 

dB _ d{Bv) ^ d^B ^ _ I d^B 
1h~ ~~2l)z^' 

In astrophysical conditions the classical resistivity is very small and the 
magnetic held behaves macroscopically as if the dihusion term in the magnetic 
induction equation would be zero. This behaviour is dictated by the magnetic 
Reynolds number, Rm = UL/rj where U and L are characteristic velocity and 
space-scale respectively. However, under certain conditions, in relatively small 
volumes, the diffusive behaviour becomes dominant and the magnetic held lines 
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reconnect. To model this in our grid, we shall assume that the physics occurs 
in two different regimes and the switch between these two regimes is dictated 
by the behaviour of the control parameter. The control parameter tells us what 
is the value of the difference between the magnetic field between one grid point 
and its neighbours. By definition this control parameter is thus local, and its 
characteristic scale I is very small. If L is the characteristic length scale of 
the simulation, then in our case Z/L is at least 0.002. If the control parame¬ 
ter exceeds a certain threshold, then what happens locally becomes worthwhile 
inspecting. Since we may assume the velocity does not change in order of mag¬ 
nitude, and since the diffusivity rj is constant (the plasma does not change), the 
ratio between the macroscopic Reynolds number and that of the local Reynolds 
number is of the same order of I jL. This can be viewed as a reason why the con¬ 
trol parameter changes the behaviour of the Reynolds number. To summarise, 
the convective behaviour is the main framework. If the control parameter be¬ 
comes critical the magnetic field evolution is given, for a brief period of time, 
by a diffusive behaviour. Once this local criticality is relaxed, control is given 
back to the convective behaviour. 

Stated more clearly, the induction equation is 

^ (4) 

1 << 1 - 


We denote the term v in 0 by S{z,t) and consider it as a stochastic 
source term. 

Equations Q are brought to dimensionless form by the transformations 
t = aT, B = bBo, z = /3Z, v = G = gBo/P"^, where a is a characteristic 
timescale, Bq is the characteristic magnetic field, /3 is the a characteristic length 
scale defined here as the distance travelled in time a by a perturbation moving 
with Alfven velocity va in a medium with a characteristic magnetic field equal 
to Bq, and wq is the initial velocity of the flow. 

With these parameters, the diffusive description becomes 

i. I, 


while the convective description becomes 





+ s{Z,T), 


avQ 



( 6 ) 


where s{Z,T) 
rameter is 


aS{z,t)/B q] the dimensionless expression for the control pa- 


1 d^b 
2 ^' 


(7) 
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3 Discretization of the induction equation 

The discretisation of the induction equation follows the general rules 


h,j+i bjj 

AT 

bi,j+i — b, 


(AZ) 


2 {bi+i,j T bi—i j j] 


AT 
9i,j = 


- = -xb^,j 

1 


- K: 




“h S 


AZ 






( 8 ) 
(9) 

2(AZ)^ ^ 

where the notation fij stands for the value of the parameter f{Z,T) evaluated 
at Z = iAZ and T = j AT. 

For the stochastic source term, at each time step j, we choose to update the 
magnetic field at a randomly chosen position kj as follows 

Si,3 = + ATe)5i^kj, ( 11 ) 

where e is some positive small number < 1. 

The diffusion description, written in update form and using Eq. (10) becomes 

h,j+i = bij - Cgij, C = 2kAT. (12) 

We will consider flows with a velocity decreasing as the spatial grid index 
increases, such that Vi+ij < Vij, and, even more, propagating flows such that 
the velocity is zero for points not yet reached by the wavefront, i.e., Ei+i j —= 
—Vij. With this assumption and written in update form, the convection is 
described as 


bi,j+i — bij [1 + VijO] + ATsij, 0 — X 


AT 

AZ' 


(13) 


When the critical threshold has been reached in a point of the grid, \gij \ > 
gcr, the critical quantity is redistributed among the neighbours of this point. 
This corresponds to the diffusive behaviour of the induction equation. 


bij+i —> bij 

which fixes ( = 2/3, and leads to a redistribution 


bi±i,j+i bi±ij + -^9i,j- 


(14) 


(15) 


The energy released by one volume during one individual relaxation event is 
the magnetic energy lost in that volume 


= 


1 

2^0 

1 

2^0 


Jv 

{Bfn - B}i„) dxdy 


f dxdy 


dz. 


(16) 
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Since the magnetic energy is not dependent on the x and y coordinates, and 
the X and y dimensions are negligible with respect to the z dimension, we take 
the result of the integral with respect to these variable to be 1. Rearranging 


and writing in dimensionless form, Eq. (161 becomes 


eR= [ {bir,-b%„)dZ, 

JdZ 


(17) 


where Er = cr (multiplied by unit area) and cr is a dimensionless quantity. 

We approximate Eq. 0 with its discretized counterpart, for one cell with 
AZ = 1, we drop the index R, and use the same notation, thus obtaining 

(18) 






This is the dimensionless energy released by one unstable cell i. The star on 
j is adopted to make the following thing clear: within one time step j of the 
simulation, the same cell, due to next neighbour interaction might become un¬ 
stable more than once. The j* is a subdivision of the simulation time step and 
it is nonzero as long as the cell is unstable. The quantity to be compared with 
observations, the energy emitted by the entire grid during one time step j is 


Cj — ^ Cij* — 

unstable nodes i 



where Nj is the number of events needed to relax the grid during time step j. 


3.1 Simulations 


The procedure is summarized below 


1. Initialization: a one dimensional grid with Nz grid points is initialized to 
hold in each cell a value for the magnetic field, bo; the initial velocity (Vb,o) 
is some multiple of the characteristic Alfven speed for the configuration 

2. Evolution: Eor each j in the interval 1,Nt time steps, the evolution of 
the system is implemented as 


choose a random number kj in the interv al l ,Nz and update the 
value of the cell in agreement with Eq. (11) (stochastic loading). 

since the upward flow with velocity V is deterministic, one can for¬ 
mally determine what cell ij the flow has reached at the current time 
step j. The numerical value in this cell is updated according to the 
first term on the left hand side of Eq. (13) (convective behaviour). 
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• at this point the mean value of the magnetic field in the grid is 
calculated; each component of the grid is then scaled with respect to 
this mean value. The critical parameter gcr (the difference between 
neighbouring cells needed to trigger a relaxation) is taken as 10% of 
the mean value. 


for each cell in the grid, i.e., i in the interval 1, Nz, gij is calculated 
with Eq. (10). If the absolute value of the control parameter is 
larger than g^r, then a diffusion behaviour is implemented in line 
with equations (141 and (15). This sweep of the entire spatial grid is 
done while cells with critical parameter can still be found. A variable 
Nj stores how many times (for one particular time step j) cells in 
the grid were relaxed. 


3. Results: Nj represents the number of events needed at each time step j 
in order to fully relax the grid. Each of these numbers is what is also 
called the avalanche size and the vector N is used to produce the event 
size distribution. 


4 Self organized criticality in one dimensional 
magnetic flows 

We discuss a set of simulations we performed, with parameters Nt = 10®, 
Nz = 500, bo = 1, e = 0.3, AT = AZ = 1, y = 6* and in which we varied the 
profile of Vi^j according to the laws Vi^j = constant, Vij = \/Vij = j~^ 
and for various values of y. For comparison purposes the results in which 
convection is ignored are also included (marked with E = 0). SOC occurs when 
cells in the system have a value for the parameter larger or equal to the critical 
parameter and when this state is spread out in the simulation grid with no 
preferred length scale. The lack of preference is seen in the power-law shape 
of the event size distribution D{N). D{N) is the number of times during the 
simulation in which a number N of events was needed to relax the grid. A fit 
of the type D{N) ~ N°‘ consistently gives a value of a close to 1 for all the 
parameter combinations shown in Table as it is expected theoretically for one 
dimensional systems. 

The event size distribution for different velocity profiles is shown in Fig. 
for y = 1. The control parameter stabilizes over the grid (e.g. in Fig. and 
the magnetic field divergence is found to be under 20%. 

The choice of parameters for the simulations set k = 1/3. This is equivalent 
to saying that rj = Ijlif?ja which is the same order of magnitude as, e.g., in 
Eq. (34) of [14] in terms of characteristic time-and-length scales. The value 
for y = VojvA sets the values of the initial velocity to that of multiples of the 
characteristic Alfven velocity. 

Samples of energy release of the grid for different velocity scalings are shown 
in Fig. right. 
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Figure 1: The evolution of the critical parameter with simulation timestep, for 
X = 5 and Vij = \/j~^- 


Table 1: Schematic display of the main characteristics a system in SOC exhibits 
(left column) and their illustration in our ID CA model. 


SOC 

ID magnetized flow 

discrete space 

ID grid 

local interaction 

one grid cell interacts with 2 other 
grid cells only at criticality (diffu¬ 
sion) 

infinitely slow external drive 

stochastic loading 

SOC occurs as a result of threshold 
dynamics 

when the critical parameter gcr is 
reached 

dissipation 

boundary dissipation 

observables 

distribution of the number of 
avalanches needed to relax one per¬ 
turbation 
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Figure 2: Left: The natural logarithm of event size distribution D{N) as a 
function of the number of events N, for different velocity scaling laws, x = 1. 
Right: The energy release for different ^velocity scaling laws, for x = 1- 
















Table 2: Parameter correlations for current simulations. The Dimensionless col¬ 
umn contains parameters which are set beforehand, and which generally char¬ 
acterise the simulation grid; the Independent parameters are those set by ob¬ 
servations; the Dependent columns contains the parameters with an analytical 
dependency with respect to the dimensionless and/or independent parameters. 


Dimensionless 

Dependent 

Resulting values 

Nt = 10^ 

P = avA 

f3 = a - 10^ m 

Nz = 500 


Afl ~ a • 4 • 10^2 J 

AT = AZ = 1 

At = aAT 

At = a ■ s 

X €{1,5,10,100} 

Az = /3AZ 

Az = a ■ 


tsira = NraAT 

tsim = a • 10® s 


Zsim = NzfdAZ 

Zsim = a • 5 • lO'^ m 

Independent 


Bo = lO^u T (obser- 
vations) 

VA = 10'^ m/s (corre¬ 
sponding to Bo) 

a 


5 Application to GRB X-Ray afterglows 

The X-Ray afterglows occur in a strongly (a half angle of 5 degrees) collimated 
magnetized outflow |35j . with initial boundary condition given by the input 
flow and mass. The flow is usually modelled as a cylinder with the lower base 
anchored on the polar region. Magnetic reconnection starts on the lower surface 
and propagates within the cylinder. We assume that this propagation occurs 
as a self organized critical (SOC) phenomenon (Table [^. According to the 
investigations of observational data by [35] and jH] , the observational signature 
of the data is of a system at criticality, namely of a one dimensional SOC. 

We apply the simulations and results presented in the previous section to 
GRBs. Estimates of numerical values of the input and output parameters can 
now be determined by analysis of the un-scaled equations. In case of GRBs 
the characteristic magnetic field may be considered to be the value inferred 
from observations, of the order of Bq = I0^°T. The dimensionless grid can be 
mapped to the physical grid by noticing that Az = /3AZ and At = aAT. The 
total simulation time is tgim = NtolAT. There is one more or less subtle issues 
here: Nt can be mapped directly to a real time when no relaxation events occur 
in the grid. If such relaxation events do occur, within one simulation time step, 
the number of real seconds increases as needed for the grid to relax. The total 
simulation length is Zsim = NzPAZ. The correlation between the parameters 
and how their quantitative interpretation changes is shown in Tables and [^ 

We propose a few parameters characterizing the GRB X-Ray afterglows to 
discriminate if whether or not the proposed model (and subsequent simulations) 
are in agreement with the observations: the time it takes the system to reach 
SOC, the energy released during the flaring, the slope of the event size distri¬ 
bution, a. Another commonly used parameter is the spectral slope (the slope 
of the power spectra) but this is shown to be in a bijective relation with a. 
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Table 3: Illustration of how astrophysical parameters (released energy, total 
duration and total length) depend on the chosen characteristic timescale (a). 


a 

Resulting values 

a = 10“"^ s 

10"« J (for = 1) 


^sim — 10 S 


Zsim = 5 • 10^ m 

a = 10"^ s 

ERr^A- 10"'' J (for CR = 1) 


tsini = 1000 S 


Zsim = 5 • 10^ m 


The time it takes the simulation to reach SOC depends on the initial loading 
of the grid. As this is generally not known, we cannot infer any quantitative 
result. 

Since the model itself does not discuss (nor can it) the direction of emis¬ 
sion, then isotropy and anisotropy of the radiation cannot be decided in this 
framework. The amount of radiation produced can be estimated based on the 
dimensionless energy unit (Table [^. 

Quantitatively, the highest dimensionless energy recorded in the simulations 
was 604000 (for constant velocity with y = 100), leading to an energy estimation 
of at least 10^® J (for a = 10“^ s). We are not in the position to offer a strict 
interpretation of this result because a clear connection between one point on 
our plots and one count in the detector does not exist. However, we can argue 
that this estimate is just a minimum: assume a sampling rate of 1 s used to 
obtain the observed light curve; further, assume that the detector records an 
integral of the emitted light curve with time as an independent variable. With 
a = 10“^, a detector set to observe our simulated light curve would observe 
at least 10^^ J. Also, notice that an increase of one order of magnitude in Bq 
would lead to an increase of three orders of magnitude in the energy estimation 
of the simulations. 

A fit of the type D{N) ^ N~°‘ for the event size distribution of all the 
simulated lightcurves are close to the value of a = 1; this is not a new result but 
merely a consistency check, as ID SOC is expected to produce slope unity in the 
event size distribution (see e.g. [T] for a derivation of this result). Simultaneously, 
this result is in agreement with results obtained from observational data by [48] . 

Qualitatively, for comparison purposes, we present sample data for GRBl 11022854 
recorded in the ng Nal detector of Fermi. We find that the bursting behaviour 
is well reproduced by a model with convective behaviour of the type H ~ 1/ t/r 
and a large ratio of initial-to-Alfven velocity (Fig|^. 

In assessing the importance of the chosen velocity profile with respect to 
reproducing observational data, there are two aspects which can be discussed 
separately: 1) reproduction of the ID event size distribution and 2) reproduction 
of observed emitted energy (light curve). As can be seen in Figurethe model 
reproduces the expected event size distribution, regardless of the velocity scaling 
law considered. However, the bursting character observed in the light curve is 
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reproduced only for the specific velocity scaling shown in Figure]^ 


6 Discussions 

There are several issues which must be stressed regarding this toy model. First, 
the model seems to show that there is a possibility that in real phenomena 
magnetic reconnection occurs in an one dimensional set-up. Second, the al¬ 
ready known problem of representation in a cellular automata: have we really 
represented the magnetic field and its dynamics? 

In real natural phenomena one dimensional magnetic reconnection is not in 
agreement with Maxwell’s, and subsequently with the MHD, equations. One 
always needs at least two independent vectors in a basis to be able to describe 
the physics of magnetic reconnection - and in fact there exist 2D MHD setups 
that do this job 

[H], m, and [20j have shown that in stochastic magnetic reconnection mag¬ 
netic field line wanderings occur in a random fashion, leading to increased re¬ 
connection velocity. They calculate the ratio of the wandering length to the 
main direction of the magnetic field and this ratio is very small. This is how 
the toy model presented in this work should be regarded. Our one dimensional 
grid obscures the fine details of the phenomenon and follow the macroscopic 
evolution of the system in agreement with some microscopic laws; and this is 
what Cellular Automata do. We place this issue in the ” sweep the microscopic 
physics under the rug” category. 

The question of representation is what can now be called an ’’old” problem. 
In a classical three dimensional setup it has already been shown that a CA 
where the represented variable is the vector potential, one can claim that both 
qualitative and quantitative data can be extracted [13 [H]. This is also valid for 
representing the magnetic field directly, with the only problem that V • H = 0 
cannot be controlled, but its validity can be checked a-posteriori, during the 
simulations. 

In short, this model does not claim to do anything more than provide a CA 
approach to a convection plus critical diffusion equation, stemming from the set 
of MHD equations written for a two-dimensional configuration. This solution, at 
least qualitatively, fits the one dimensional Self Organized Criticality signature 
inferred from observations and theory. 


7 Conclusions 

Following the suggestion by [48] that X-ray flares in GRBs originate from mag¬ 
netic reconnection - driven events, we have investigated the possibility of forma¬ 
tion of Self-Organized Criticality in an one dimensional magnetized fluid flow. 
As mentioned in [ISJj such a work ”... could not only help to understand the 
central engines of GRBs, but also help to study applications of solar magnetic- 
reconnection theories.” In the present paper we have adopted a simplified theo- 
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retical approach for the study of magnetic reconnection, which is based on the 
one-dimensional form of the magnetic induction equation, with the background 
flow explicitly included. In the induction equation we have considered the term 
—vdB/dz as a stochastic source term, and we have constructed an equivalent 
cellular automaton model to describe the evolution of the magnetized plasma 
flow, as well as its energy emission. 

Our analysis, based on computing the event size distribution and the fact 
that a grid averaged value of the critical parameter stabilises to a constant 
value during the simulations, concludes that SOC appears in a one-dimensional 
magnetized setup with background flow. This result was not unexpected, as 
the onset of magnetic reconnection is threshold dependent and many observa¬ 
tional and theoretical efforts, together with simulations have established that 
observational data can be explained by a self-organizing (dynamical) spread of 
reconnection events. 

Even more, SOC seems to be an underlying principle of Nature. One novel 
pillar in this line of thought has been set by the observational analysis performed 
by [IH], and sustained by previous efforts to explain the diversity of GRBs 
through one single principle. This is why the possibility of the application of 
the present model to the GRBs afterglows was also briefly considered. Both 
qualitative and quantitative results show that the same grid, equipped with the 
same evolution law is able to produce a large diversity of light curves (e.g., 
Fig.l), while conserving a power-law shape for the event size distribution. 
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Figure 3: Upper figure: Light curve for GRB111022854 recorded with the ng Nal 
detector (see [SI mi). Lower figure: The energy release for a velocity scaling 
of V ^ l-y/r and y = 10. 
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